Illuminating Dark Matter at the ILC 



Herbert Dreiner* 

Physikalisches Institut and Bethe Center for Theoretical Physics, University of Bonn, Bonn, Germany 

Moritz Huck^ and Michael Kramer^ 
Institute for Theoretical Particle Physics and Cosmology, RWTH Aachen University, Aachen, Germany 

Daniel Schmeier§ and Jamie Tattersall 1 

Physikalisches Institut and Bethe Center for Theoretical Physics, University of Bonn, Bonn, Germany 

(Dated: November 9th, 2012) 

The Wimp (weakly interacting massive particle) paradigm for dark matter is currently being 
probed via many different experiments. Direct detection, indirect detection and collider searches 
are all hoping to catch a glimpse of these elusive particles. Here, we examine the potential of the Ilc 
(International Linear Collider) to shed light on the origin of dark matter. By using an effective field 
theory approach we are also able to compare the reach of the Ilc with that of the other searches. 
We find that for low mass dark matter (< 10 GeV), the Ilc offers a unique opportunity to search 
for Wimps beyond any other experiment. In addition, if dark matter happens to only couple to 
leptons or via a spin dependent interaction, the Ilc can give an unrivalled window to these models. 
We improve on previous Ilc studies by constructing a comprehensive list of effective theories that 
allows us to move beyond the non-relativistic approximation. 



I. INTRODUCTION 

Weakly interacting massive particles (WlMPs) are one 
of the leading candidates to solve the dark matter puzzle 
[1] . Primarily this is due to the fact that a neutral parti- 
cle that interacts with roughly the strength of the weak 
force, naturally gives the correct relic abundance. In ad- 
dition many theoretical models predict that the masses of 
these states should exist around the scale of electroweak 
symmetry breaking, e.g. Supersymmctry (Susy) [2, 3], 
Universal Extra Dimensions (Ued) [4], Little Higgs [5] 
etc. 

Currently, this Wimp paradigm is being actively ex- 
plored in a number of different ways. Perhaps the most 
well known are the direct detection searches that aim 
to observe interactions between the dark matter and an 
atomic nucleus [6]. As these are extremely low rate ex- 
periments, the detectors are typically placed deep under- 
ground to reduce background. The annihilation of dark 
matter into Standard Model particles in high density re- 
gions of our universe offers another potential method to 
see a signal e.g. [7]. 

In particle colliders here on Earth the same interac- 
tions may be probed in the production of dark matter. 
Unfortunately, the fact that WlMPs are neutral and only 
weakly interacting means that they cannot be detected 
directly in these experiments. Therefore collider based 
searches must rely on particles produced in combination 
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with the dark matter candidates. If dark matter is pro- 
duced directly, one possibility is to use initial state radia- 
tion (Isr), such as gluon jets, or photons, that will recoil 
against the WlMPs. 

This idea was first explored in a model independent ap- 
proach for the International Linear Collider (Ilc) using 
mono-photons in a non-relativistic approximation [8, 9]. 
Later, detailed detector studies have been performed to 
understand the full capabilities of the Ilc for such a 
signature [10-14]. Furthermore the same signature has 
been considered in the case of Susy [15, 16]. At the 
Lhc (Large Hadron Collider) and Tevatron similar sig- 
nals have also been studied but with a mono-jet signal 
[17-25]. All of these papers used the idea of parameter- 
ising the dark matter interactions in the form of effective 
operators. This has the advantage that the bounds can 
be compared with those coming from direct detection 
and also that a non-relativistic approximation is not re- 
quired to compare with the relic density measurement. 
These methods have now been used by the Lhc experi- 
ments to set bounds on different effective operators that 
are competitive with other methods [26, 27]. In addition, 
Lep (Large Electron-Positron Collider) data has been re- 
interpreted to determine corresponding constraints [28]. 

In this paper we take the effective field theory approach 
to dark matter and apply this to an Ilc search [29-31]. 
To apply the effective field theory in a consistent way we 
assume that the dark matter particles can only interact 
with the Standard Model fields via a heavy mediator. 
The mediator is always assumed to be too heavy to be 
produced directly at the Ilc and thus can be integrated 
out. For our model choices we consider the possibility 
that the dark matter candidate could be a scalar, a Dirac 
(or Majorana) fermion or a vector particle. The same 
choices are taken for the heavy mediator and all combi- 
nations are considered. The collider phenomenology can 
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vary significantly, depending on whether the mediator is 
exchanged in the s- or t-channel and consequently we ex- 
amine both. In addition, we also study the different ways 
in which the mediator can couple to both the dark matter 
and Standard Model particles. We note that using the 
effective field theory approach allows us to move away 
from the non-relativistic approximation that had previ- 
ously been used in Ilc studies. This can be especially 
important if the dark matter candidate happens to be 
light. 1 

For all models we compare the reach of the Ilc with 
the bounds derived from direct and indirect detection. 
We also calculate the couplings expected to lead to the 
correct relic density and see whether the Ilc can probe 
these regions of parameter space. We also note that an 
Ilc search is complementary to that at the Lhc thanks 
to the different initial state. 

The paper is laid out as follows. We begin in Sec. II 
by explaining how we derive the effective field theories 
for the dark matter interactions and we explicitly give 
the Lagrangian for both the full and effective theory. We 
also describe the benchmark models that we use through- 
out the study. In Sec. Ill we describe the various astro- 
physical constraints on our effective theories. We begin 
with the calculation of the relic density abundance before 
moving on to explain the bounds from direct and indirect 
detection. 

Section IV describes in detail the potential search for 
dark matter at the Ilc. Here we explain the calculation 
of the signal rate and the dominant backgrounds that 
were considered. In addition we detail how the Ilc de- 
tectors are modeled to account for relevant experimental 
effects. We find that the polarisation of incoming beams 
is particularly important for many models of dark mat- 
ter to discriminate the signal and background. We also 
investigate the advantage of a doubling of the Ilc energy 
to y/s = 1 TcV. 

In Sec. V we present the results of the paper. We begin 
by examining the potential bounds of the Ilc on the ef- 
fective coupling of the dark matter model at the collider. 
Afterwards, we combine these results with those from 
direct and indirect detection to understand for which 
models and mass ranges the Ilc presents a unique op- 
portunity to discover dark matter. Finally in Sec. VI we 
conclude and summarise the main results of our work. 



II. MODELS 

A. General Motivation 

The idea of parametrising the interaction of a dark mat- 
ter particle with Standard Model particles by using effec- 
tive operators is not new, see for example Refs. [18, 28, 



1 The mass determination of a light ncutralino dark matter candi- 
date at the Ilc has been discussed in Ref. [32]. 



30, 31, 33, 34]. Many authors construct a list of effective 
4-particle-interactions with Lorentz-invariant combina- 
tions of 7^, and spinor-/ vector-indices up to mass 
dimension 5 or 6. In many cases there is no explanation 
how those operators may arise in an underlying funda- 
mental theory. That makes it difficult to judge how ex- 
haustive the lists of operators are, whether interference 
between different operators should be taken into account 
and how the effective model is connected to realistic fun- 
damental theories and their couplings. 

We follow the effective approach introduced in [31] by 
starting from different fundamental theories with given 
renormalisable interactions between Standard Model 
fermions and the hypothesized dark matter particles that 
are mediated by a very massive particle. From these 
theories we deduce effective 4-particle-vertices for ener- 
gies significantly smaller than the mass of the mediator. 
Working with these effective operators, one can deduce 
information about the effective coupling and propagate 
this information to the parameters of the corresponding 
underlying fundamental theory. The effective approach 
allows us to reduce the dimensionality of the parameter 
space and more easily compare the different experimental 
searches. 



B. Deriving Effective Lagrangians 

We start with a list of fundamental Lagrangians 
taken from [31]. However we do not perform a non- 
relativistic approximation, since we are interested in the 
phenomenology of this Lagrangian at a high energy ex- 
periment and therefore the results for our effective oper- 
ators differ. We also use a different method to evaluate 
the effective vertices, motivated in Ref. [35], which uses 
the path integral formalism. 

We give one explicit example for the derivation of the 
effective operators and only mention specific peculiarities 
for the other cases, which are apart from that calculated 
similarly. Let ip be a Standard Model fermion and \ a 
complex scalar field representing the dark matter candi- 
date. For our example, we assume the mediator to be a 
real scalar field, </>, with mass Mq (we will keep this no- 
tation for the mediator mass throughout). The relevant 
terms in the UV completed Lagrangian are then given 

by 

^uv - \ [d^(xf - \m^ 2 (x) - g xX \x) X {x)4>{x) 
- Hx) (g s + ig P i b ) i>(x)(p{x) , (1) 
= - l -<j>{x)U x ^{x) - \m& 2 (x) - F(x)<t>(x) . (2) 

where the function F{x) is given by, 

F(x) = . 9x x f (x)x(x) + i>{x) (g s + ig pl b ) $(x) . (3) 

We have not included the kinetic terms for \, ip, as they 
are not relevant for the computation of the effective La- 
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grangian. In this particular example, g s , g p are dimcn- 
sionless couplings and g x is a dimension one parameter 
but these definitions can change depending upon the pre- 
cise model studied and we shall use this notation through- 
out. We have included the kinetic term for <p, the heavy 
mediator field. After integrating out <j>, we obtain the 
effective Lagrangian, 



off 



1 

2MH 



f 2 d j^xh4>(g s + ig P j 5 )^- (4) 



Cases with different spin for the dark matter or the me- 
diator particle are evaluated similarly. We only want to 
give some special remarks: 

• For spin-1/2 mediators, the Dirac propagator has 
only one power of Mq in the denominator, 



1 



Mr, 



1 



M 2 n 



(5) 



for these. We also choose two possible values for Mq 
to represent different suppression scales of the respec- 
tive second order terms. Models with real fields that are 
trivially connected to the corresponding complex cases 
by multiplicative prefactors are not taken into account 
separately. We also omit models with left-handed cou- 
plings that are related to the respective right-coupled 
cases. Information on these can easily be extracted from 
the related models by rescaling the corresponding result 
accordingly. 



III. ASTROPHYSICAL CONSTRAINTS 

Any model which aims to describe dark matter, for ex- 
ample through a Wimp, has to agree with present data. 
It has to give the correct relic abundance, and must be 
consistent with the bounds from direct and indirect de- 
tection searches [37-39]. 



We therefore get two effective vertices after expand- 
ing the Lagrangian up to order 1/Mq. 

• Some effective operators give derivatives on the 
Standard Model fermion fields. These are not neg- 
ligible, since they only vanish if the Dirac equation 
i0ip = mip can be used and the fermion mass m is 
small. This is not the case for e.g. heavy quark con- 
tributions in the annihilation sector and processes 
with off-shell fcrmions. 

• We use the same list of effective operators for the 
cases of real scalar (x = x 1 ), real vector (xn = xV) 
or Majorana fermion [36] dark matter fields. How- 
ever, we would like to mention that for consistency 
we do not introduce additional factors of 1 /2 in the 
couplings as is often done in the case of real fields. 

The full list of models with their respective fundamen- 
tal and effective Lagrangians is given in Table I. Note 
that all Lagrangians are hermitian by construction. 



C. Benchmark Models 

The effective operators described above have multi- 
ple independent parameters to describe the effective cou- 
pling, for example g x , gi, 9r and Mq in the scalar dark 
matter, vector mediator (SV) case or g s ,g p and Mq in 
the fermion dark matter, scalar mediator (FS) case in Ta- 
ble I. Considering the full range of parameters would lead 
to a plethora of scenarios, well beyond the scope of this 
paper. Thus we restrict our analysis to specific bench- 
mark models (see Table II) with constraints on the in- 
dividual couplings such that only one overall multiplica- 
tive factor remains. The effective coupling constant G for 
each model is then defined as G = gigj/M^. For models 
with fermionic mediators, the leading term has only a 
1/Mq dependence, which is why we define G = gigjfMQ 



A. The Relic Abundance 

We first consider the best measurement of the relic 
abundance from Wmap-7 [37] , 



DM.2 



0.1099 ± 0.0056. 



(6) 



We employ the solution of the model dependent Boltz- 
mann equation obtained in [30] , 



1.04-10 9 GeV" 



x f 



m P iy/g4xf){a + 3b/ x f ) 

(7a) 



Here mp\ is the Planck mass. Xf = M x /Tf is the inverse 
freeze-out temperature, Tf, rescaled by the Wimp mass, 
M x . It is implicitly given by the equation, 



x f = In 



c(c + 2) 




g m p \M x {a + 6b/xj) 



• (7b) 



g*(xf) denotes the relativistic degrees of freedom in equi- 
librium at freeze-out and is given in Ref. [40] . a and b are 
the first two coefficients of the non-relativistic expansion 
of the thermally averaged annihilation cross section, 



(av) w a + bv 2 + 0{v 4 ), 



(8) 



where v is the relative velocity of the colliding particles. 
Here the center-of-mass energy squared is approximated 
by [33,34], 



4M^+My+3/4JWy. 



(9) 



g are the internal degrees of freedom of the Wimp, c is an 
order unity parameter which is determined numerically 
in the solution of the Boltzmann equation and we set this 
parameter to 0.5. 
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DM Med. Diagram 



V 



V 



V 



F V 



F tS 



F tV 



V V 



-.zuv 

-J&eB 



9x ,,t 
Ml 



(9s - gp^i'x'x + ttxH (g 2 s + g 2 P - ^gsg P i 5 ) 7"9 M (i>x) 




Mn 




/ 9x(.^d^x ~ xd^)Z» + H>Y(9iPl + grPn^Z^ 



ffx j. 
Ml 



^7" (giPL + g r P R ) V {fd^ - 4>d^) 




x {gs + g P i 5 ) x<t> + ^ {gs + g P i'°) ^<\> 
-j-pX {gs + ig P i 5 ) xi> {gs + ig P i 5 ) i> 




WiaiPz. + 9rPR)i>Z^ + X1"(9iPl + grPR)xZ„ 

^2 ^7" (9\Pl + grPn) ip x7m [giPL + grPo) x 

X (gs + g P r') i><P + 4< (g s + g P j 5 ) x4> 
jpip {gs - 3p7 5 ) xx (g s + 3p7 5 ) ^ 

WiaiPt + 9rPR)xZ» + xY(giP L + BrPaWz? 
^7 M (giPL + grPn) xx7m (giPL + grPn) ip 




x^xA + Tp(gs + ig P -y 5 )ip4> 
||x M x„^ (gs + ig P i°) i> 



m^giPL + grP R )x» + i>^(giPL + grP R )vxl 
1 gigr^Y^P xlx P + -^xl^Yi^Y (sf-Pt + 9tPr) 3 M (i>xp) + 




Wx [Z^xldx^ + ZpXvdxT + xlXvdZ^] + ^(g,P L + g r P R )^ 



Wx J„ M 

Ml 



^7 M (giPL + grPo) ip \x V 9xU ~ X^^Xm" + d " (xlxi* ~ xlxA] 



TABLE I. List of interaction vertices for S(calar), F(ermion) and V(ector) dark matter, x, before and after integrating out 
the heavy mediator scalar field cf>, spinor field rj or vector field with mass Mn- denotes the Standard Model fermion. 
dX^ v = d^X v — d" V . tS and tV denote cases where the mediator is exchanged in the t-channel. 



Operators 


Definition 




Name 


SS, VS, FS, 


g P 


= 




scalar 


FtS, FtSr: 


gs 


= 




pseudoscalar 


SF, SFr: 


g P 


= 0, Mn = 


1 TeV 


scalar _low 




9p 


= Q,Mn = 


10 TeV 


scalar -high 




9s 


= 0, Mn = 


1 TeV 


pseudoscalar Jow 




9s 


= 0, Mn = 


10 TeV 


pseudoscalar -high 


SV, FV, FtV, 


gi 


= gr 




vector 


FtVr, VV: 


91 


= -g r 




axialvector 




gi 


= 




right-handed 


VF, VFr: 


gi 


= g r ,Mn = 


1 TeV 


vectorJow 




gi 


= —g r ,Mn 


= 10 TeV 


vector -high 




91 


= g r ,Mn = 


1 TeV 


axialvector -low 




91 


= —g r , Mn 


= 10 TeV 


axialvector_high 


FVr : 


gi 


= 




right-handed 



TABLE II. Benchmark models with specific values for the 
coupling constants shown in Table I. 



Instead of testing all the models presented in Table I, 
we shall focus on a few exemplary cases. First, the relic 
density depends on the possible Standard Model parti- 
cles, /, the WlMPs can annihilate into XX ~~ * f f ■ We 
shall consider two cases for the set of particles /: (i) 
all leptons, (ii) all SM fermions. Second, two variants 
of couplings are tested. In one scenario all SM parti- 
cles couple via the mediator to the Wimp with the same 
strength; this is called universal coupling. In the other 
they have a coupling proportional to their mass, which 
we call Yukawa-like coupling. In the cases where we 
have the same effective operator our results agree with 
Refs. [33, 34], up to the normalisation (see Appendix A). 

In order to set constraints, we must determine the total 
relic density, which is the sum of the relic density of the 
particle and the anti-particle (if the latter exists) . This 
means the relic density for a complex particle-pair is two 
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times the density of a real particle. If we consider the 
Wmap result as an upper bound on the relic density, i.e. 
allowing for other dark matter, then this corresponds to 
a lower bound on the effective coupling of the Wimp to 
the SM particles. If we require our Wimp to be the only 
dark matter, we shall also obtain an upper bound on the 
effective coupling. 

The strict interpretation that our model only contains 
a heavy mediator and a single Wimp ensures that there 
are no resonances or co-annihilations. However we also 
note that in many full theories that contain dark matter, 
a 'co-annihilation' regime can exist that can significantly 
alter the relic density in the universe. Whilst the co- 
annihilation mechanism cannot be incorporated into the 
strict definition of our model, it may actually have no ob- 
servable effect on the collider based phenomenology. An 
example of such a feature could be stau co-annihilation in 
Susy that would not change the Ilc production process 
of the lightest supersymmetric particle. Another example 
is that a more complicated model may contain resonant 
annihilations. Both of these examples can significantly 
weaken the relic abundance bounds. 



B. Direct Detection 

We shall also impose bounds on our operators from the 
direct detection searches for Wimp dark matter. The ex- 
periments are designed to measure the recoil energy from 
the scattering between a (dark matter halo) Wimp and 
the target nucleus. The interactions are difficult to detect 
since the energy deposited is quite small, 1 to 100 keV, 
[1]. These experiments give an upper limit for the cross 
section between the dark matter and the nucleus of the 
target. One drawback is that in the cases where the 
Wimp does not couple to quarks, the coupling can only 
occur through loop diagrams. 

The direct detection experiments give a much stronger 
bound on spin independent (SI) interactions than on spin 
dependent (SD). The reason is that in the SI case the 
interaction with all nucleons add coherently which en- 
hances the corresponding cross section by the atomic 
number squared. However, the spins of the nucleons can- 
cel if they are paired. Thus SD interactions are only 
enhanced for very special nuclei. 

The SI interactions are scalar or vector interactions in 
the s-channel, the axialvector and tensor interactions in 
the s-channel give a SD interaction. Note that due to the 
low kinetic energy of the Wimps the cross section should 
be computed in the non-relativistic limit. In that case 
the pseudoscalar interaction, ip^ip, vanishes. 

The i/u-channel diagrams are cast into a sum of s- 
channcl diagrams via the Fierz identities. From this only 
the SI parts are employed, since any SD contribution is 
negligibly small. Tensor interactions occur only via the 
Fierz identities, since we do not consider fundamental 
tensor interactions. However, since Fierz identities will 
always give at least one SI contribution, tensor terms can 



be dropped. 

For the SI interactions we shall consider the limits set 
by the XenonIOO experiment [38]. These are the most 
recent and set the strictest limits over a broad parameter 
range. For the SD interactions we consider the XenonIO 
data [41] since XenonIOO gives no statement on SD in- 
teractions. The smaller data set along with the physical 
reasons mentioned above lead to a bound that is <~ 10 6 
times weaker than for the SI interactions. The calcula- 
tions for the WlMP-nucleus cross sections follow Ref. [31] 
and for identical models we find the same results. See 
Appendix B for the complete list of cross sections. 

C. Indirect Detection 

We also consider the indirect detection searches for 
dark matter. These are much more model dependent, 
as the dark matter is seen via an agent, for example neu- 
trinos, which could also be produced via other means. 
Specifically we shall consider the Pamela experiment 
[39] which measured an excess of positrons. These could 
potentially originate from dark matter annihilation. To 
implement this we need to compute the propagation of 
the produced positrons and electrons from the source to 
the earth. This is described by the diffusion-loss equa- 
tion [42], 

^ - V[tf(x,£)VV] - ^[b(EW = q(x,E). (10) 

Here ip(x, E) = dn e+ /dE is the positron density per en- 
ergy. K(x,E) is the diffusion coefficient which describes 
the interaction with the galactic magnetic field. b{E) de- 
notes the energy loss due to synchrotron emission and 
inverse Compton scattering. q(x, E) is the source term 
due to dark matter annihilation. We note that convec- 
tion and re-acceleration terms are ignored as these do 
not apply to positrons [43]. 

We use the conventional formalism [44, 45] to derive 
a solution of Eq. (10). It is also possible to use the so- 
called extended formalism that takes the corrections from 
sources in the free propagation zone into account as well 
as those from the diffusion zone. However, this increases 
the runtime of the calculation considerably while only 
giving a small correction that is less than the measure- 
ment error. To perform the numerical comparison we use 
the cored isothermal dark matter density profile [46] and 
the galactic propagation model M2 [44]. 

The above choices result in the following positron flux, 

$ e+ (E) = ^(r Q ,z Q ,E), (11) 
i,(r, z,E)= T -^ J*™* de s f(es)I(r, z, e, e s ), (12) 

it \ Y^V^ t ( OLiT \ ■ n7r ( z + L ) 
I(r, z,e,es) = 2_ J 2^ J ^~j^l sm 2L 

i n 

x exp(-u iin (t-ts))Ri,n, (13) 



() 



(14) "0 



Here rg, R, Kg, L are parameters which describe 
the M2 propagation model. They are set to the stan- 
dard choices [44, 45] t e = 10 16 s, R = 20 kpc as well 
as to the M2 propagation model L = 1 kpc, Kq — 
0.00595 kpc 2 /Myr, 5 = 0.55. /(e) is the energy distri- 
bution of the positrons from the annihilation and is gen- 
erated with Pythia8 [47]. are the coefficients of the 
Bessel- Fourier expansion of R(r, z), 



R(r,z) 
p(r, z) 



r){crv) 



Po(— ) 7 
r 



1 + (r Q /r s ) a 
1 + (r /r) a 



(0-j)/a 



(15) 

(16) 



Here (an) is the thermally averaged annihilation cross 
section. We include all possible final states, not just those 
resulting in positrons. Furthermore ry = 1/2 for real par- 
ticles and 1/4 for complex particles. r = 8.5 kpc is 
the distance of the solar system from the galactic center. 
Pq — 0.3 GeV/cm 3 is the local dark matter density and 
a = /? = 2, 7 = 0, rs — 5 kpc are chosen according 
to the cored isothermal dark matter density distribution 
[44, 45]. 

Pamela measures the ratio $ e +/(<l> e - + < E > e +), where 
the fluxes, $ e ± , contain the flux from dark matter an- 
nihilation and from any astrophysical background. The 
background we take is [42] , 



e bg 



dE 



0.16e 



-l.i 



+ 



1 + lie 09 + 3.2e 215 
0.7e 07 



1 + HOe 1 - 5 + 600e 2 - 9 + 580e 



4.2 



e+bg 



4.5e l 



0.7 



dE 1 + 650e 2 - 3 + 1500e 4 - 2 



GeV^cixrW" 1 , 
(17a) 

GeVW-VW -1 , 

(17b) 



e = E/GeV. 
The quantity we compare to Pamela is, 



4>, 



ei-x 



e+bg 



4> 



e+bg 



e-bg 



(18) 



and we note that $ e + x = ®e- x - 

We find an upper bound on the annihilation cross sec- 
tion by assuming that all of the excess comes from dark 
matter. However, it is possible that other background 
sources contribute and thus we also allow models that 
produce a flux smaller than the one seen. 

We also note that for dark matter masses above 
~1 TeV, the Fermi-Lat [48] experiment may provide 
competitive bounds from inverse Compton scattering 
[49, 50]. However, since we are only interested in models 
that can be probed at the Ilc we ignore them here. 
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FIG. 1. Diagrams for radiative pair production of dark mat- 
ter. Terms in which the heavy mediator can emit a photon 
are neglected. 



The IceCube collaboration also sets limits on heavier 
dark matter masses via annihilations into neutrino final 
states [51, 52]. In addition these bounds may be com- 
petitive for spin dependent interactions but we do not 
consider the limits in this study. 



IV. DARK MATTER SEARCH AT THE ILC 

A. Radiative Production of Dark Matter 

For the Ilc search, we look at the process e + e~ — » XX7 
with a hard photon being the only detected particle in the 
final state, Fig. 1. We determine the polarized differential 
cross section for this process with respect to the relative 
photon energy x = 2£ , 7 /y / s and its polar angle 9 by 
integrating over the full phase space of the final state dark 
matter particles. The results for this calculation are given 
in Table III, with further explanation of the abbreviations 
used given in Appendix C 1. Previous Ilc studies, e.g. [8, 
13], have used the Weizsacker- Williams approximation 
for soft photons. This formula relates the differential 
photon cross section to the total pair production cross 



section e + e~ — 
s — » s = s(l 
function F x g, 

da [e + e~ 



XX with a reduced center of mass energy 
— x) and multiplied by the kinematical 



XX7] 



da; dcos0^ 



Ft 



a e T e -> xx\ 



(19) 



Due to the soft collinear approximation used, we expect 
that the above equation will perform poorly for large an- 
gle and high px photons. We compare the analytical 
result to this approximation to test the reliability. In Ta- 
ble III we put terms in bold, which are purely caused by 
our analytical treatment. The corrections are either of 
the form of an additional kinematical factor V x g, mostly 
appearing in models with vector mediators, or completely 
new terms that typically appear in t-channel interac- 
tions. Since lim a; _j.o V x g — 1 and lmx^o^i) = 0, the 
WW-approximation is in agreement with our full result 
for small energies. In Fig. 2 we show the respective pho- 
ton energy distributions for different models in both the 
WW-approximation and the full analytical treatment. 



7 



Model 

da; dcosfl 



32^M£ 
SF ^ 



2 PS 4 

G 2 _ P C S + 12M 2 [(ff s + 3p) 4 C H + (g s - g P ) 4 C L ] + A SF 



sv |^^?^+ 5 2 c.] 5 2 



FV [G i+r s/3 2 + 3 (g, + g r f M 2 ] [gfC L + g 2 C R ] 



1 n 

FVr (s, -grf [gfC L + g 2 C R ] 



FtS 4§^tG? +p - Af 2 ) + A FtS ] 



FtSr J^_G 2 S+P [3(1 - 2M 2 )C P + K e 2(s - 4M 2 )Cy] 
FtV [6G? r C s (s - 2M 2 ) + (§- M 2 )VM9?Cl + g 4 r C R )] 



FtVr 



[l2G 2 r C s (s - 2M 2 ) + {s — 4M 2 )V x e(9tC L + g 4 r C R )] 



VF SsJImiM* [ 40G ^^( 7M x - + s 2 ) + ^ (gfC L + g A r C R ) (40M X 6 - 22M A J + 56M 2 s 2 + 3s 3 ) + A VP ] 

VFr SsJJmiM* [ 60G 'rC*( 12M x ~ *M 2 s + s 2 ) + ^ (sfCi + 5r 4 C«) (320M X 6 - 104 4 J + 32M 2 s 2 + s 3 ) + A VFl ] 

W 3j^ ^ + ggj ffx(-^x + 20M 2 J + 

TABLE III. Analytical differential cross sections for the process e + e~ — ¥ \X1 m the various effective models. Terms in bold do 
not appear in the Weizsacker- Williams approach and are given in Appendix C 1 where we also define all used abbreviations. 
Models with a suffix 'r' correspond to the case of real particles. Cross sections for SSr, FSr and VSr are twice as large as in 
the complex case while SV and VV vanish completely for real particles. 



The curves behave quite congruently with differences vis- 
ible in the high energy sector. Since most of the signal 
events lie in the low energy part, the approximation gives 
accurate results for counting experiments. A shape de- 
pendent analysis would need to use the analytical result 
to estimate the correct threshold behaviour for high en- 
ergies. Our subsequent analysis is performed with the 
full analytical cross section. 

When we restrict the various couplings in our model 
according to the benchmark scenarios, Table II, most of 
the cross sections simplify and have only one polarisation 
dependent term Cj. To determine the polarisation lead- 
ing to the best signal to background ratio, we only need 
to consider cases with different Cj. We therefore classify 
our models as follows: 

Scalar-like : cr po i = Csa u „ P oi, (20) 
Vector-like : er po i = Cycr unpo i, 



Right-like : cr po i = C R <r unpo \, 
Left-like : cr pol = C L a unpol . 

Models with t-channel mediators usually have multi- 
ple terms with different polarisation behaviour and do 
not fall into one of the basic polarisation classes given in 
Eq. (20). We choose the following polarisation settings 
for those: 

• Models with fermionic mediators are classified ac- 
cording to their leading term, which is always 
scalar-like. 

• All other models have both scalar-like and vector- 
like parts of about the same size. We analyse them 
in a vector-like scenario that naturally leads to a 
better background suppression. 




FIG. 2. Comparison of tree level photon energy distributions in the WW-approximation and the analytical solution for 
M x = 50 GeV, | cos6> 7 ! max = 0.98 and = 1 TeV. (a) SV, (b) FtS. 



B. Standard Model Background for Monophotons 

We consider the two leading dominant Standard Model 
background contributions after selection, determined 
with a full Ild (International Linear Detector concept) 
detector simulation [13, 53]. All numbers here and in the 
following paragraphs refer to the nominal Ilc center of 
mass energy of 500 GeV [54] . We also consider the case of 
an increased energy of 1 TeV and mention the differences 
later. 

• Neutrinos from e + e~ — > form a polarisa- 
tion dependent background. The leading contribu- 
tion comes from t-channel W -exchange, which only 
couples to left-chiral leptons. Additional smaller 
contributions come from s-channel Z-diagrams 
with both left- and right-chiral couplings. We also 
consider the case of one additional undetected pho- 
ton, which contributes with a size of roughly 10 %. 

• Bhabha scattering of leptons with an additional 
hard photon, e + e~ — > e + e~7 has a large cross sec- 
tion but a very small selection efficiency, since both 
final state leptons must be undetected. It has been 
determined to give a contribution of the same or- 
der of magnitude as the neutrino background, after 
application of all selection criteria. It is mostly po- 
larisation independent [13, 53]. 



C. Data Modeling 

To evade the use of a full detector simulation, we 
build on the results of Refs. [13, 53]. For the signal 
and monophoton neutrino background, we generate the 
events by ourselves with the given phase space criteria. 
We then apply the Ild estimates for the energy resolu- 
tion as well as the reconstruction and selection efficien- 
cies 2 and compare the final energy distributions. For the 
diphoton neutrino and Bhabha background, we model 
the final distributions directly from the given results per- 
formed with a full detector simulation [13, 53]. 

For the generation of signal and monophoton neutrino 
events we use CalcHEP [55]. We produce signal events for 
all benchmark scenarios with dark matter masses ranging 
from 1 GeV to 240 GeV. To avoid collinear and infrared 
divergences, we limit phase space in the event generation 
to E 7 e [8 GeV, 250 GeV] and cos6» 7 e [-0.995,0.995]. 
Initial State Radiation (Isr) and beamstrahlung signif- 
icantly change the width and position of the neutrino 
Z°-resonance, Fig. 3a), and are taken into account. We 
set the accessible parameters in CalcHEP according to 
the Ild Letter of Intent [56] to 645.7 nm for the bunch 
size, 0.3 mm for the bunch length and a total number of 
particles per bunch of 2 • 10 10 . 

The finite resolution of the detector components and 
the use of selection criteria to reduce beam-induced back- 
ground are taken into account by applying the following 
steps to both signal and background data. First we shift 



Other background sources contribute with less than 1 % 
compared to the neutrino background and are therefore 
omitted. 



2 From here on, the expression 'efficiency' abbreviates 'reconstruc- 
tion and selection efficiencies'. 
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FIG. 3. Photon energy distribution before and after application of beam effects (Isr + beamstrahlung) and detector effects 
(resolution + efficiency) for a) unpolarised neutrino background and b) unpolarised FS scalar signal with M x = 150 GeV. 
Distributions are normalised to 10 6 tree level events. 



□Signal 
□ Background 

vvy 

vvyy 

e + e"y 




250 
E y in GeV 

FIG. 4. Photon energy distributions of the most dominant 
background contributions (stacked) compared to an example 
signal (FS Scalar, M x = 150 GeV) with a total cross section of 
100 fb. All spectra are taken after selection for an unpolarised 
initial state. 



p- 


/P+ 




1/1/77 


e + e 7 


+0.8/ 
-0.8/ 


0/0 
-0.3 
+ 0.3 


2257 (2240) 
493 (438) 
5104 (5116) 


226 (228) 
49 (43) 
510 (523) 


1218 (1229) 
1218 (1204) 
1218 (1227) 



FIG. 5. Total number of events in the different background 
sources after application of all selection criteria. The numbers 
are given for an integrated luminosity of 1 fb -1 in different 
polarisation settings. Numbers in brackets are taken from 
Ref. [53] which employed a proper detector simulation. 



the photon energy, given in GeV, according to a Gaus- 
sian distribution by taking into account the estimated 
resolution of the Ild detector components [56], 



16.6 % 



E 1 ^/S 7 in GeV 



1.1 %. 



(21) 



Afterwards we further limit the phase space to re- 
duce background processes in the Z° resonance peak at 
242 GeV and additional collinear photons from Isr, 



£ 7 e [10 GeV, 220 GeV] , 
cos6> 7 € [-0.98,0.98] . 



(22) 



The additional angular cut ensures a good photon recon- 
struction within the detector. Finally a random elimina- 
tion of events is used to simulate the efficiency factor for 
reconstruction and selection determined in Ref. [53] . The 
efficiency consists of an energy dependent part t\ and a 
constant part 62 that are applied successively, 



£2 



97.22 % 
96.8 %. 



(E~ in GeV) ■ 0.1336 %, (23) 



Fig. 3 shows how these settings affect the signal and back- 
ground spectrum and Fig. 4 shows a stacked histogram of 
the dominant background processes along with a example 
dark matter signal. 



D. Analysis 

We are interested in determining the lower bound on 
the effective coupling constants that the Ilc can find for 
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P-/P+ 


iV B 


A stat 
^50 


A stat 
^500 


A sys A sys 


A tot 

^50P 


A tot 

50P 


X tot 
^500P 


X tot 
^500P 


0/0 


184998 


















+0.8/+0.3 
+0.8/+0.6 


97568 
102365 


312 
320 


99 
101 


312 
385 


125 
154 


441 
500 


336 
355 


327 
398 


159 
184 


+0.8/-0.3 
+0.8/-0.6 


87974 
83177 


297 
288 


94 
91 


169 
104 


68 
42 


341 
307 


304 
291 


193 
138 


116 
100 


-0.8/+0.3 
-0.8/+0.6 


341597 
404970 


584 
637 


185 
201 


351 
501 


140 
200 


682 
811 


601 

668 


396 
546 


232 
284 


-0.8/-0.3 
-0.8/-0.6 


212851 
148478 


461 

385 


156 
122 


233 
337 


93 
135 


517 
512 


471 
408 


275 
359 


173 
182 



TABLE IV. Total amount of background events, Nb, with statistical error, A s at , systematic error, A sys , and the total error, 
A ° . The subscripts 50 and 500 denote the integrated luminosity in inverse femtobarn. In case of a ten times larger luminosity, 
one will get ten times as many events in all channels; to better compare to the error of the low luminosity case, we give 
Ajfcfo* = Af a 7\/l0. The polarisation uncertainties are set to 0.25 % (P) and 0.1 % (P). 



each individual model under the assumption that no sig- 
nal events are measured. We perform a counting exper- 
iment by using the TRolke [57] statistical test. We 
determine the total number of background events along 
with its statistical and systematic fluctuation ANb and 
exclude coupling constants which would lead to a larger 
number of signal events than the 90 % confidence interval 
of the background-only assumption. 



E. Systematic Uncertainties 

Systematic uncertainties play an important role in de- 
termining the total error on the background, ANb, and 
for estimating the bounds on the effective couplings. 
There are two dominant contributions, motivated in 
Ref. [53] which we now discuss. 

The experimental efficiency given in Eq. (23) will be 
determined at the real experiment by measuring the Z°- 
rcsonance peak, which is theoretically known to a very 
good accuracy. Systematic uncertainties on that value 
are given by the finite statistics of this measurement and 
further broadening of the peak by unknown beam effects. 
These errors can be extrapolated down to the dark mat- 
ter signal region at small photon energies and, since the 
same efficiency factor is used for signal and background, 
is highly correlated between those two. This global un- 
certainty will therefore approximately cancel in the de- 
termination of the maximum coupling G e ff. 

Cancellation will not take place for model dependent 
effects however. This is due to the fact that the sig- 
nal energy distribution depends on the unknown mass of 
the dark matter particle and the underlying interaction 
model. Therefore, the correct function e{E 1 ) for the sig- 
nal will be different from the used neutrino background 
efficiency given in Eq. (23). Since we do not know the 
model a priori, we use the same value for both and intro- 
duce an error on the determination of the signal events, 
Ns- Compared to Ref. [53], we use a conservative value 
of Ae = 2 %. 

Since the neutrino spectrum depends on the incoming 



lepton's polarisation P ± , any fluctuation within those 
parameters will give additional systematic uncertainties 
on the number of expected background events. One can 
not use the information from measuring the Z° resonance 
in this case to infer information in the low energy sig- 
nal range because of the polarisation dependence of the 
shape itself. Given the assumed accuracy of at least 
AP/P = 0.25 % [56] with a possible improvement to 
0.1 % at the Ilc, we can derive the corresponding error 
on the polarised number of background events. As an 
example we show the left handed background, 

N pol = (l+P+)(l-P-)N unpoh 

AN pol = ^[P-(l + P+f + [P+(l-P-)} 2 ^ AUpoi- 

(24) 

From the numbers in Table 5, we assume an identical 
polarisation dependence for uvj and 1/1/77 events and 
no dependence for the Bhabha background. 



I A type P~/P+ 


N s 


»"50P 


r 50P 


?*500P 


r 500P 


Scalar +0.8/ + 0.3 


620 


1.41 


1.85 


1.90 


3.90 


+0.8/ + 0.6 


740 


1.48 


2.08 


1.86 


4.02 


Vector +0.8/ - 0.3 


620 


1.82 


2.04 


3.21 


5.34 


+0.8/ -0.6 


740 


2.41 


2.54 


5.36 


7.40 


Left -0.8/ + 0.3 


1170 


1.72 


1.95 


2.95 


5.04 


-0.8/ + 0.6 


1440 


1.78 


2.16 


2.64 


5.07 


Right +0.8/ - 0.3 


1170 


3.43 


3.85 


6.06 


10.09 


+0.8/ -0.6 


1440 


4.69 4.95 


10.43 


14.4 



TABLE V. Determination of the best ratio r = N S /AN B 
with AJVb given by the different total errors determined in 
Table IV. Ns describes the number of polarised signal events 
for the different classes described in Sec. II with a common 
reference value of 500 unpolarised events for an integrated 
luminosity of 50 fb _1 . We only show the polarisation signs 
with the largest ratios. We mark the numbers which lead to 
the best signal to background ratio in bold. 
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P-/P+ 


iV B 




^500 


51 


*? 


A tot 

^50P 


A tot 

^sop 


A tot 

^500P 


X tot 
^500P 


0/0 


162437 


















+0.8/+0.3 
+0.8/+0.6 


54649 
62791 


234 
251 


74 
79 


380 
469 


152 
188 


446 
531 


279 
314 


387 
476 


169 
203 


+0.8/-0.3 
+0.8/-0.6 


38365 
30223 


196 
174 


62 

55 


201 
125 


82 
50 


281 
214 


212 
181 


210 
137 


102 
74 


-0.8/+0.3 
-0.8/+0.6 


357173 
435979 


598 
660 


189 
209 


428 
612 


171 

245 


735 
900 


622 
704 


468 
647 


255 
322 


-0.8/-0.3 
-0.8/-0.6 


199561 
120755 


447 
348 


141 
110 


284 
411 


114 

165 


530 
538 


461 
385 


317 
425 


181 
198 



TABLE VI. Total amount of background events (Nb) and different error sources (see Table IV) for ^Js = 1 TeV. 



F. Polarisation Settings 

Polarisation can be used to significantly increase the 
number of signal events according to Eq. (20) but also 
increases the systematical contribution to the total back- 
ground error, AN-g via Eq. (24). We are interested in 
the settings for each individual model that leads to the 
largest N S /AN B ratio allowing for the strictest bounds 
on G c ff- In Table IV we give the total number of back- 
ground events in different polarisation settings P~ = 
±0.8 and P+ = ±0.3/±0.6 that are feasible at the Ilc 
[54] . We give the statistical fluctuation for integrated lu- 
minosities of 50 fb _1 as well as for 500 fb _1 . Since the 
latter will give ten times as much events in all channels, 
we reduce the statistical error accordingly to give a value 
comparable to the small luminosity case. We also give 
the systematic error that is dominated by the polarisa- 
tion uncertainty for two estimates of the polarisation er- 
ror SP/P = 0.25 % and 0.1 % [58]. Finally we give the 
total errors adding all combinations of individual errors 
in quadrature. 

On the signal side, we look at the different classes de- 
rived in Sec. II with respect to their polarisation depen- 
dence. For comparison, we use a common reference value 
of 500 unpolarised events for an integrated luminosity of 
50 fb _1 and derive the corresponding number of events 
for polarised input. 

We look for the maximum ratio r = N$ / ANb and the 
results for the best settings are displayed in Table V. In 
most cases the largest possible polarisation for the incom- 
ing leptons enhances the result. For high statistics and 
a non-reduced polarisation error, the systematic uncer- 
tainty from increased polarisation may outweigh the gain 



P-/P+ 






e + e 


0/0 


2677 


268 


304 


+0.8/ -0.3 


421 


42 


304 


-0.8/ + 0.3 


6217 


622 


304 



TABLE VII. Simulated and modeled number of events in the 
different background sources after application of all selection 
criteria for ^/i = 1 TeV. The numbers are calculated for 
an integrated luminosity of 1 fb _1 in different polarisation 
settings. 



in the number of signal events though. In those cases, 
which appear only in scalar- and left-coupling models, 
less polarised beams lead to better results. 

G. Increasing a/s to 1 TeV 

We also consider the possibility of a doubled center 
of mass energy. This changes the previous analysis as 
follows: 

• We generate events in a larger photon energy range 
Ery <s [8 GeV, 500 GeV] and reduce it to the inter- 
val [10 GeV, 450 GeV] after performing the energy 
resolution shift AE/E. This again reduces back- 
ground events from the Z°-resonance, which now 
is positioned at 496 GeV. 

• Dark matter signal processes can now be produced 
with masses up to 490 GeV. 

• We use our previously modeled distribution for the 
Bhabha background and rescale it by a factor of 
1/4, taking into account that the full cross section 
for that process is approximately proportional to 

1/8. 

• We use, as a rough approximation, the same Isr- 
and beamstrahlung parameters in CalcHEP, effi- 
ciency factors and systematic error estimates. 

Tables VI- VIII summarise again the number of back- 
ground events per background scenario, the individual 



Model 


P-/P+ N s 


?*50P 


r 50P 


^BOOP 


r 500P 


Scalar 


+0.8/+0.3 620 


1.39 


2.22 


1.60 


3.7 




+0.8/+0.6 740 


1.39 


2.36 


1.55 


3.65 


Vector 


+0.8/-0.3 620 


2.21 


2.92 


2.95 


6.08 




+0.8/-0.6 740 


3.46 


4.09 


5.40 


10.00 


Left 


-0.8/+0.3 1170 


1.59 


1.88 


2.50 


4.59 




-0.8/+0.6 1440 


1.60 2.05 


2.23 


4.47 


Right 


+0.8/-0.3 1170 


4.16 


5.52 


5.57 


11.47 




+0.8/-0.6 1440 


6.73 


7.96 


10.51 


19.46 



TABLE VIII. Determination of the best ratio r = N S /AN B 
(see Table V) for a/5 = 1 TeV. 
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error sources and the determination of the best polarisa- 
tion setting for the increased center of mass energy. In 
contrast to the Bhabha cross section that falls mainly 
according to a oc 1/s, the neutrino background gets sig- 
nificant contributions from t-channel W ± s, which give 
s/m^y -terms in the evaluation of the total cross sec- 
tion. The left-handed neutrino contribution therefore 
gets enhanced whereas the Bhabha background becomes 
less dominant in some polarisation channels. This leads 
to a larger relative polarisation error and therefore a 
larger impact on the size of the background fluctuation. 
In the end, vector- and right-coupling models receive 
stronger enhancement for polarised input than in the 
\A = 500 GeV case, whereas the other models suffer from 
the larger impact of polarisation on the total error and 
prefer smaller polarisation. 



V. RESULTS 

We begin by presenting the reach at the Ilc in terms 
of the effective coupling constant in Sec. VA. We then 
compare these potential bounds with the couplings pre- 
dicted by the cosmological relic density and the bounds 
coming from direct and indirect detection experiments. 
Of course we would also like to discover a dark matter 
at the Ilc and the bounds provide an estimate of the 
potential sensitivity of the collider. 

A. ILC Bounds 

We determine the 90 % exclusion bound for the ef- 
fective coupling constant in each benchmark model for 
the best case scenario. The integrated luminosity is 
set to 500 fb -1 and the systematic polarisation error to 
AP/P — 0.1 %. For each benchmark model we choose 
the polarisation setting that leads to the best signal to 
background ratio for the corresponding polarisation be- 
haviour according to Tables V and VIII. Results for dif- 
ferent polarisation settings can be found by res calin g the 
bound on the coupling according to G' = Gsjr'/r with 
r denoting the ratio Nq/ANb given in Table VIII. We 
choose to present all of the results for an Ilc with a cen- 
ter of mass energy of 1 TeV due to the increased range 
of dark matter masses that this option can probe. In ad- 
dition, smaller effective couplings can be probed, mainly 
due to the falling Bhabha background. 

In Fig. 6 we show the derived bounds on the coupling 
constants for an Ilc center of mass energy of 1 TeV. The 
hashed area denotes the region that either violates the 
tree level approach with a too large dimensionless cou- 
pling constant g 2 > iir, or by having a too small mediator 
mass Mq < 1 TeV, for the effective approach to be valid. 
Note that the leading order in models with fermionic 
mediators has a different mass dimension and therefore 
gives a different definition for the effective coupling con- 
stant Goff- If a model has no separate 'pseudoscalar' or 



'axialvector' results, it is identical to the corresponding 
'scalar'/ 'vector' line due to identical cross section for- 
mulas. For masses away from the threshold, the Ilc is 
able to exclude coupling constants down to the order of 
10" 7 GcW 2 or 10" 4 GeV" 1 , depending on the mass di- 
mension. This corresponds to a total cross section (for 
the given phase space criteria) of about 0.3 fb. Excep- 
tions however arise for models with vector dark matter 
that tend to have very strong exclusion limits for small 
masses. This is caused by the 1/M 4 dependence in the 
photon cross section, which leads to divergences for very 
small vector boson masses. It has been shown [59] that 
only spontaneously broken gauge theories can lead to 
models with massive vector particles that are not diver- 
gent. Therefore, our initial fundamental model cannot 
be the full theory for all energies. In our effective ap- 
proach, we restrict the energy to a maximum and in that 
case one can still receive perturbative valid results for 
mass ranges that do not violate unitary bounds. How- 
ever, the perturbativcly allowed mass range cannot be 
given in this model independent approach, since such an 
analysis needs more information about the size of the in- 
dividual couplings and the relation between the mass of 
the mediator and the dark matter mass itself. In sum- 
mary, a more detailed fundamental theory is needed to 
evaluate the breakdown of perturbation theory in this 
scenario. 

We note that in models with fermionic operators, the 
sub-leading order has a negligible effect, as can be seen 
from the nearly identical lines for fermionic mediators 
with different masses. 

B. Combined Results 

The combined maximum exclusion limits for spin in- 
dependent DM-proton interaction at Pamela, Wmap 
and the Ilc are shown in Figs. 7-9. We choose a subset 
of models that couple to all Standard Model fermions 
and give an overview of the bounds that we can ex- 
pect. Other models behave similarly and are therefore 
not shown again separately. We can give the following 
statements about the comparison of the Ilc exclusion 
bound with the current XENON limits: 

• We have sensitivity to spin independent proton 
cross sections for, as an example, the FV Vector 
model down to 10" 42 cm~ 2 or cquivalcntly 10" 3 fb, 
which is an improvement of about four orders of 
magnitude compared to current Lep [28] and two 
orders of magnitude compared to current Tevatron 
[18] and Cms [26] results. 

• An increased center of mass energy can lead to 
stronger bounds by up to one order of magnitude. 
It also allows a larger dark matter mass range to 
be probed. 

• Ilc bounds get significantly weakened if the inter- 
action is Yukawa-like. At the Ilc the mediator 
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FIG. 6. 90 % exclusion limits on the effective couplings accessible at the Ilc with y/s — 1 TeV. We only give effectively allowed 
regions for models with dimensionless fundamental couplings g. 
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must couple to electrons, which have a suppressed 
Yukawa coupling. The production cross section is 
thus small, leading to weaker bounds. 

• Models with scalar mediators give weaker bounds 
than models with vector interactions. For fermionic 
dark matter we observe a difference of about two or- 
ders of magnitude, which is in agreement with pre- 
viously mentioned results from e.g. Lep. For scalar 
and vector dark matter the difference is mass- 
dependent and can increase to up to six orders of 
magnitude, which is due to the different mass di- 
mension of the couplings. 

• The Wmap bounds are for many effective models 
very constraining, Figs. 7-10. However, we would 
like to point out that these can be highly depen- 
dent on the full theory whilst not affecting the Ilc 
or direct detection phenomenology. For example, 
annihilation can occur via some resonance or as in 
some Susy models, co-annihilation with staus or 
stops. 

In Fig. 10 we show some models which allow for lepton 
couplings only. In that case, dark matter can only in- 
teract with protons via photons through a fermion loop, 
cf. Appendix. B 4. The loop factor significantly lowers 
the cross section and therefore increases the bound in the 
case of vector coupled models. Other models allow quark 
couplings only at the two-loop level or theoretically com- 



pletely forbid them [28]. In all cases, the Ilc would give 
the strongest exclusion bounds for dark matter lepton 
couplings. For models with fermionic mediators there is 
an extra subtlety when comparing the bounds. In partic- 
ular the exclusion limit at the Ilc is mainly given by the 
leading term in the operator expansion, which is scalar 
like. Loop couplings can only happen for vector currents, 
which in the case of a fermionic mediator is only given 
by the sub-leading order and has an additional factor of 
1/Mq. In that case, when translating any exclusion lim- 
its into bounds on the WiMP-proton cross section, we 
need to know the exact mass of the mediator. We show 
this in Fig. 10 for the two different chosen suppression 
scales 'Low' (Mq =1 TcV) and 'High' (M n =10 TcV), 
Table II. 

In Fig. 11 we show the exclusion limits for the spin- 
dependent interaction. In our case, only the model with 
fermionic dark matter, a vector mediator and an axial- 
vector coupling leads to such an interaction. In that case, 
we compare with data from the previous Xenon exper- 
iment (XenonIO), since no results for the XenonIOO 
phase were available when this study was completed. 
Since in this scenario dark matter only couples to a 
single nucleon on average because of the natural spin 
anti-alignment in nuclei, the Xenon bounds are not 
coherently enhanced by the atomic number and there- 
fore strongly lose sensitivity. The Ilc would also give 
strongest exclusion bounds over the whole accessible 
mass range here. 
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FIG. 7. Combined 90 % exclusion limits on the spin independent dark matter proton cross section from Ilc, Pamela and 
Wmap for a selection of scalar dark matter models. 
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FIG. 8. Combined 90 % exclusion limits on the spin independent dark matter proton cross section from Ilc, Pamela and 
Wmap for a selection of fermionic dark matter models. 
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FIG. 9. Combined limits on the spin independent dark matter proton cross section from Ilc, Pamela and Wmap for a selection 
of vector dark matter models. 
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VI. CONCLUSIONS 

In this paper we considered a broad range of effective 
models for dark matter and investigated the possibility 
that these models could be explored at the Ilc. The 
models considered the possibility that dark matter was 
a new scalar, fermion or vector particle and would be 
produced at the Ilc via a new, heavy intermediate state, 
the mediator particle. For the mediator we also consid- 
ered spins 0, 1/2 and 1. We obtained the corresponding 
effective theories by integrating out the mediator field. 

To be able to compare the reach of the Ilc with the 
other experimental searches, certain assumptions have to 
be made on how the mediator and dark matter couples 
to the Standard Model particles. We assume in all mod- 
els that interactions only occur with the Standard Model 
fcrmions but the relative strength to different particles is 
varied. In the simplest variant we choose that the cou- 
pling is equal between all the Standard Model states. An- 
other choice is that the interaction scales with the mass 
of the interacting Standard Model fermion, a 'Yukawa- 
like' interaction. The last choice we make is the most 
optimistic for Ilc phenomenology with only the Stan- 
dard Model leptons interacting with the heavy mediator. 
Since the produced dark matter particles will be invisi- 
ble to the Ilc detectors, we require a radiated photon to 
be emitted from the initial state that will recoil against 
missing momentum. This topology provides a distinctive 
signal with which to discover dark matter. For the Ilc 
study, we included the dominant backgrounds and most 
important detector effects. In addition we considered the 
possibility of using polarised initial states to reduce back- 
grounds and improve the signal strength. 

The effective theories that we consider provide an ef- 
ficient way to compare the reach of the Ilc with other 



methods to discover dark matter. Firstly, we consider 
the dark matter annihilation cross section required for 
the relic density observed by Wmap. We also look at 
the direct detection bounds at Xenon by calculating the 
dark matter-nucleon scattering cross section. In addi- 
tion, we include bounds from dark matter annihilation 
to positrons from the Pamela experiment. 

In terms of the effective dark matter model, we 
found that the Ilc should be able to probe couplings 
If/ 7 GcV~ 2 , or 10" 4 GcV- 1 depending on the mass di- 
mension of the theory. In models that contain vector 
dark matter, the Ilc may be able to probe even weaker 
couplings in the case of low dark matter mass. 

To compare with astrophysical bounds, we found that 
the Ilc reach is strongly dependent on the exact dark 
matter model. If we assume that dark matter is rela- 
tively heavy (> 100 GeV) and interacts with a Standard 
Model particle in proportion to its mass, then the Ilc is 
uncompetitive. However, in the case that dark matter is 
relatively light (< 10 GeV) then the bounds from the Ilc 
arc competitive with astrophysical bounds in many mod- 
els. In addition, if dark matter happens to only interact 
with the Standard Model leptons then the Ilc offers a 
unique possibility to discover dark matter. For this rea- 
son, an Ilc search is complementary to those done at the 
Lhc thanks to the different initial state. 
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Appendix A: Cross Sections for Annihilation 

We give the full cross sections for annihilation of a 
pair of dark matter particles with mass M x into a pair 
of Standard Model fermions with mass m,f. To find the 
expansion coefficients in av rts a + bv 2 , we perform the 
non-relativistic approximation s w 4M 2 +M 2 v 2 -|-|M 2 v 4 
[30]. Note that in order to find the correct result for 
the v 2 term in av, it is necessary to expand up to order 
v 4 because of the appearance of yfs in the cross section 
formulae. 

The total cross section is then given as the sum of the 
cross sections over all allowed final state fermions. This 
set is restricted both by kinematics (m/ < M x ) and by 
the assumed model. The latter also determines whether 
the coupling Gf is universal or particle-dependent. 

We define the mass ratio £ = rrif/M x and the velocities 
of both particles fix = \A ~ 4m| / s to compactify the 
following expressions. 

Some of our effective operators have been analysed be- 
fore, for example [33, 34], and we agree with the respec- 
tive results for the annihilation cross sections. 
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Appendix B: Cross Sections for Direct Detection 

We now give results for the dark matter-nucleon scat- 
tering cross section at zero momentum transfer, <j°, for all 
defined benchmark models. In a universal scenario, the 
effective coupling is independent of the quark (G q = G), 
whereas it grows proportionally to the quark mass in a 
Yukawa- like model (G q = G m q /m e ). We use the follow- 
ing definitions: 
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with the numerical values for fP and A? listed in [60, 61]: 



fP = 0.020 ± 0.004, 


(B5) 


f p = 0.026 ± 0.005, 


(B6) 


ff = 0.118 ±0.062, 


(B7) 


A p u = -0.427 ±0.013, 


(B8) 


A^ = 0.842 ±0.012, 


(B9) 


Af = -0.085 ±0.018. 


(B10) 



The cross sections can be evaluated in a nonrelativistic 
approximation for the Wimp and by using the quark pro- 
ton form factors listed above. Sec e.g. [31]. If a model 
is not listed, its scattering cross section equals zero, e.g. 
for pseudoscalar interactions that always vanish in a non- 
relativistic model. Again, we agree with the respective 
results in [33, 34] for comparable operators. 

Cross sections for real final state particles can easily 
be derived from the following list by setting the vector 
form factors b p and b p to zero and rescaling f p and d p by 
a factor of 2. 
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Furthermore we define the reduced mass of the Wimp 
proton system, 
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4. Photon Loop 

If the Wimp only couples to leptons, the WlMP-proton 
interaction can only happen at the loop level. In that 
case, a low energy photon that couples to a virtual lep- 
ton pair interacts with the whole proton. This only hap- 
pens for models with s-channel vector bilinears V>7^, i.e. 
models which include either a, b p , or a, b p , term in the 
low energy tree level cross section. Results can therefore 
be derived as follows, 
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where the reduced cross section has to be understood 
as the tree level cross section given above after setting 
bp, b p = 1 and f p , d p = 0. This ensures that we only take 
the vector interaction parts. If the tree level cross section 
includes a b p term, the loop factor is given as, 



I 

For bp terms, it reads, 
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In both cases, the loop function can be evaluated as, 
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We follow the conservative assumption of a maximum 
scattering angle to find q 2 = — A[i 2 v 2 with /j, describ- 
ing the reduced mass of the WlMP nucleus system and 
v = 500 km/s being the typical escape velocity of a Wimp 



in a dark matter halo. Because of the new g-dependence 
of the cross section and the fact that the photon only cou- 
ples to the protons inside the nucleus, the official Xenon 
results have to be rescaled according to, 
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where q — q{M^ = Mp) uses the reduced mass n of the 
Wimp proton system instead. This weakens the cross 
section limits by about a factor of 10. 



Appendix C: Differential Cross Section for 

e+e~ -> XXI 

1. Abbreviations 

We use the following abbreviations for the final cross 
section list in Table III: 

Polarisation prefactors: 

C s = l + P + P~, C v = l-P+P-, (CI) 
C L = (1 - P")(l + P + ), C R = (1 + P~)(l - P+). 



Terms with combined couplings: 

Gx±y = 9x ± ffy! Gxy = gxgy- 
Relativistic velocities: 



/3: 



4M 2 



4M x 
s{\ - x) 



Kinematical functions: 

a (x - l) 2 + I 



F x g 
V x g 



it x sin o 

x 2 cos(26») + (3a: 



4 ((x- 1) 2 + 1) 



(C2) 

(C3) 

(C4) 
(C5) 



We show terms that arise in the analytical evaluation of 
the differential photon cross section in e + e~ — > \X1 but 
not in the Weizsacker-Williams approximation in (C6)- 
(C10). They all vanish in the soft-photon limit x — >• 0. 



■ I - / = {1 t^- [(9s + 9 P ) 4 C R + (g s g P ) 4 C L ] 



4M 2 1 - x 



AsFr 



a s 



8tt M 2 1 - x 



[{9s+g P ) A C R + {g s -gp) A C L ] 



(C6) 
(C7) 



24 



iFtS 



(1 - V xe ) 



C s (s - AM 2 X ) + — C S (2M 2 + s) 



(C8) 



A VF = 20G 2 r C s (l - V^^is 2 + 4M 2 s - 8M 4 ) + 



(gtC L +g?C R ) 



M 2 



1 x 4 sin 2 (2(9) s (3s 2 + 26M 2 s - 32M 4 ) 
32 (x-l) 2 ((x-l) 2 + l) 



+ 6 * -s(s 2 + 7M 2 s - 24M*) - 1(1 - 14 e )(21 S 3 + 282M 2 s 2 - 1144M 4 s + 160M^) 

{{X L) + 1) 4 

| 3 (1-148) ^,2 2SM 2 S j 16M 4 ) I 1 ^— Y^1 ^C7^ iocm2; , oni\/r4\ , ( 1_ ^)^/'52 , .,,,2- , ; i ; 1 



2 {l-x 
{gtC L +g 4 C R ) 



Ml 



x ' 4 (l-x) 2 
x 4 sin 2 (26») 



s(7^ - 126M 2 s + 32M*) + ) _ Tg s(« + 2M x* + 6M x) 

(1 x) 



(C9) 



32 (x-l) 2 ((x-l) 2 + l) 



-s(s 2 + 32M 2 S - 24M 4 ) + 2- -s(s 2 + 12M 2 s + 56M 4 ) 

1)2 + 1) X X ((X — l)^ + 1) * * 



1(1 - 14e)(7s 3 + 144M 2 s 2 - 168M 4 s + 1280M^) + 1^1 ^s(s 2 - 48M 2 s + 56M 4 ) 

^ A. 2(1 3/ ) 

+ t ^"^ ^ 2 - 272M 2 s + 104M 4 ) +2 ^ - V *°} s(S 2 + 2M 2 s + 6M 4 ) 



4 (l-x) s 



(1 - x) 3 
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